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of  Pfaffian  systems.  In  addition,  we  give  new  conditions  for  converting  Pfaffian 
systems  into  Goursat  form  after  prolongation,  and  for  linearizing  control  systems 
using  dynamic  extension. 

Several  examples  of  mobile  robots  are  examined,  and  it  is  shown  that  for  some 
kinematic  arrangements,  a  prolongation  corresponding  to  a  dynamic  state  feedback 
is  needed  to  transform  the  corresponding  Pfaffian  system  into  Goursat  form. 


This  research  was  supported  in  part  by  the  NSF  under  grant  IRI-9014490  and  by  the  ARO  under 
grant  FT-DAAL03.  D.  Tilbury  would  also  like  to  acknowledge  an  AT&T  Ph.D.  Fellowship  for 
financial  support  of  this  work. 


1 


2 


TILBURY  AND  SASTRY 


1.  Introduction 

There  has  been  a  great  deal  of  interest  in  the  use  of  exterior  differential  systems 
for  analyzing  nonlinear  control  systems.  We  bring  together  here  some  of  the  results 
which  have  been  recently  published  as  well  as  add  some  of  our  own  contributions  to 
this  area.  We  show  that  all  of  the  main  results  in  exact  linearization  of  nonlinear 
systems  can  be  restated  in  terms  of  exterior  differential  systems,  and  in  addition,  we 
present  a  new  set  of  sufficient  conditions  for  linearization  by  dynamic  extension. 

In  this  paper,  we  will  only  make  use  of  a  special  type  of  exterior  differential  system 
called  a  Pfaffian  system,  and  all  of  our  definitions  and  results  are  specific  to  this  case. 
Loosely  speaking,  a  Pfaffian  system  is  represented  by  a  codistribution  of  one-forms 
edefined-on  the  state  space,  which  we  will  assume  to  be  a  connected  manifold.  These**#** 
one-forms  may  represent  constraints  on  the  system  velocities,  as  in  the  case  of  mobile 
robots  where  the  wheels  roll  without  slipping.  A  control  system  with  state  x  €  Rn, 
input  u  6  Rm  of  the  form  x  =  f(x,  u)  can  also  be  written  as  a  Pfaffian  system,  with 
the  constraints  defined  as  the  one-forms  o'  =  dxi  —  /‘(x,  u)dt  on  Rn+m+1.  Although 
special  care  must  be  taken  for  control  systems  to  treat  time  differently  from  the  state 
and  input  variables,  we  will  show  that  the  main  results  on  transforming  Pfaffian 
systems  into  normal  forms  will  carry  over  into  the  special  case  of  control  systems. 

For  mechanical  systems  with  linear  velocity  constraints,  such  as  mobile  robots  with 
wheels  that  roll  without  slipping,  exterior  differential  systems  are  the  most  appro¬ 
priate  method  for  analysis.  Using  methods  from  exterior  differential  systems,  we 
can,  under  certain  conditions,  transform  the  Pfaffian  system  defined  by  the  rolling 
constraints  into  a  normal  form  and,  in  the  new  coordinates,  easily  find  solution  tra¬ 
jectories  for  the  system. 

The  outline  of  this  paper  is  as  follows.  First,  we  present  some  background  mate¬ 
rial  on  exterior  algebra  and  Pfaffian  systems.  We  then  give  the  definitions  for  the 
Goursat  normal  form  and  extended  Goursat  normal  form.  After  stating  necessary 
and  sufficient  conditions  for  converting  Pfaffian  systems  into  these  normal  forms,  we 
show  that  several  examples  of  mobile  robot  systems  satisfy  these  conditions.  We  then 
consider  systems  which  do  not  satisfy  the  conditions  for  this  conversion,  and  present 
the  concept  of  prolongation  of  a  Pfaffian  system.  We  give  sufficient  conditions  for 
converting  a  Pfaffian  system  to  extended  Goursat  normal  form  using  a  specific  type  of 
prolongation,  and  we  show  an  example  of  a  mobile  robot  system  which  satisfies  these 
conditions.  Finally,  we  turn  our  attention  to  control  systems,  expressed  as  Pfaffian 
systems.  Since  the  Brunovsky  form  is  a  special  case  of  Goursat  normal  form,  all  of 
the  results  for  converting  Pfaffian  systems  to  Goursat  form  can  be  specialized  to  give 
conditions  for  converting  control  systems  to  Brunovsky  linear  form.  We  show  that 
the  special  type  of  prolongation  we  proposed  in  this  paper  is  the  dual  of  dynamic 
extension,  or  adding  integrators  to  some  of  the  input  channels.  Thus  our  theorem  for 
converting  Pfaffian  systems  to  extended  Goursat  normal  form  using  prolongations  can 
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be  specialized  to  give  sufficient  conditions  for  linearizing  control  systems  by  dynamic 
extension. 


2.  Exterior  Algebra  and  Pfaffian  systems 

We  give  a  brief  overview  of  some  of  the  definitions  of  exterior  algebra  and  Pfaffian 
systems  that  we  will  use  in  this  paper.  The  theory  of  exterior  differential  systems 
is  powerful  enough  to  analyze  solutions  of  systems  of  partial  differential  equations; 
however,  we  will  restrict  ourselves  in  this  paper  to  systems  of  first-order  ordinary 
differential  equations.  We  encourage  the  interested  reader  to  consult  the  monograph 
by  Bryant  et.  al.  [1],  from  which  most  of  this  introductory  material  was  taken,  for 
more  details. 

A  real  vector  space  V  or  its  dual  (covector)  space  V *  generates  an  exterior  algebra 
with  the  exterior  or  wedge  product  defined  by 

q  A  =  —fl  A  a 
a  A  a  =  0 

aa  A  (6/?  +  C7)  =  ( ab)a  A  (3  +  ( ac)a  A  7 

for  all  <*,/?€  V(V),  a,6,c  €  ft.  The  wedge  product  of  two  vectors  is  called  a  two- 
vector.  We  define  A2(V)  as  the  space  of  two-vectors.  We  can  similarly  build  up  higher 
vectors  and  define  Ak(V)  as  the  space  of  all  k- vectors.  For  completeness,  we  define 
A°(V)  =  R  and  =  V.  The  dimension  of  A*(V)  is  (”).  From  the  axioms,  it 

follows  that  Ak(V)  is  empty  for  k  >  n. 

The  exterior  algebra  over  V  is  a  graded  algebra, 

A(V)  =  A°(V)  ©  AX(V)  ©  •  •  •  ©  An(V) 

Any  element  A  €  A(V)  can  be  written  uniquely  as 

A  =  Ao  +  Ai  +  •  •  •  +  An 

where  A ;  €  A'(V)  for  t  =  1,. . .  ,n. 

Now,  consider  a  differentiable  manifold  M  of  dimension  n  and  its  cotangent  bundle 
T*M .  We  construct  the  bundle  A(TaM)  whose  fibers  are  the  exterior  algebra  of  T*M, 
that  is: 

a  (t;m)  =  a  °(r;)  ©  a \t;)  ©  A2(r;)  ©  •  •  •  ©Ancr;) 

The  bundle  A(T'M)  has  Ap(TmM)  as  sub-bundles.  A  section  of  the  bundle 

AP(T*A/)  =  U  A p(T;M) 

*€M 

over  M  is  called  an  exterior  differential  form  of  degree  p  or  simply  a  p-form. 
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For  local  coordinates  on  M  denoted  by  x  =  (a?i, . . .  ,x„)  ,  a  local  basis  for  TXM  is: 

{—  -2-j 

We  denote  its  dual  basis  on  T*M  as 

{^li  •  •  •  j 

defined  by 


In  terms  of  these  local  coordinates,  a  p-form  v  can  be  written  as 
WU-ip(X)  dx'i  A  • '  •  A  dxip  «!  <  *2  <;••_<  ip 

where  the  coefficient  functions  are  smooth  functions  on  M. 

We  will  use  the  notation  Qp(Af)  to  mean  the  module  (over  the  ring  of  smooth 
functions)  of  all  smooth  sections  of  AP(T*M),  and  Q(M)  =  ©  QP(M)  as  the  module 
of  forms  on  M. 

We  begin  by  considering  a  codistribution  I  on  M ,  spanned  by  s  one-forms,  that  is 

I  =  span{o;1,...  ,q*} 
where  a’  is  in  f ^(M)  for  *  =  1, . . .  ,s. 

Definition  1.  Pfaffian  Systems. 

On  a  manifold  of  dimension  n,  a  Pfaffian  system  is  the  smallest  ideal  X  C  fl(Af), 
generated  by  a  codistribution  I  of  one-forms  spanned  by  {o1, . . .  ,  a'}  which  is  closed 
under  wedge  products. 

Any  element  of  X  can  be  written  in  the  form: 

a  =  Y^0'Aai 

i=l 

where  O'  is  any  element  in  Q(M).  Throughout  the  course  of  this  paper  we  will 
deliberately  confuse  the  notation  and  refer  to  the  codistribution  I  as  the  Pfaffian 
system. 

The  dimension  of  a  Pfaffian  system  is  defined  to  be  s,  the  number  of  independent 
one-forms  which  generate  it.  Any  n  —  s  linearly  independent  one-forms  which  are 
independent  of  I  form  a  complement  to  I.  The  codimension  of  I  is  n  —  s. 

An  integral  curve  for  a  Pfaffian  system  is  a  curve  c(t)  :  (— e,  e)  — ►  M  which  satisfies 
the  constraints,  that  is,  c*(o')  =  0  for  all  o’  €  I.  Here  the  notation  c*(o’)  is  taken 
to  mean  (a',  ^)). 

A  (local)  independence  condition  for  a  Pfaffian  system  is  a  one-form  r  which  does 
not  vanish  on  integral  curves,  that  is  c*(r)  ^  0.  We  add  the  additional  condition  that 
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r  be  integrable,  so  that  our  Pfaffian  systems  will  correspond  to  systems  of  first-order 
ordinary  differential  equations. 

We  will  need  the  notion  of  congruence  modulo  a  Pfaffian  system.  For  I  =  {a1, ...  ,  a*} 
we  say  that »/  =  £  mod  /  if 

=  {  +  Aa< 
t=i 

for  some  forms  O'  in  f l(M). 

Definition  2.  Exterior  Derivative. 

The  exterior  derivative  is  defined  as  the  unique  map  d  :  fl k(M)  — ►  which 

satisfies  the  following  properties: 

(1)  For  /  € 

df  =  -~-dx  i  +  •  •  •  4-  •7r—dxn, 
ux\  cten 

relative  to  a  local  coordinate  chart,  or  the  usual  gradient. 

(2)  For  a€fT,£€fP, 

d{a  A  /?)  =  da  A  /?  +  (— l)ra  A  d/3. 

(3)  <P  =  0. 

Definition  3.  Derived  Flag. 

Given  a  Pfaffian  system  1  =  {a1, . . .  ,  a*},  the  derived  flag  is  defined  to  be  the  nested 
chain  of  codistributions  given  by  =  I  and 

j(*+i)  =  {a,  €  /(*) ;  dw  =  0  mod  /<*>} 

The  construction  is  assumed  to  terminate  at  some  N,  when  fl**)  =  The 

derived  flag  is  then  defined  to  be  the  sequence  of  nested  codistributions, 

/  =  J(0)  D  J(1)  D  •  •  O  T 

We  will  assume  that  the  dimension  of  jW  is  well-defined  for  all  k. 

Remark  1  ( Maximally  Nonholonomic).  The  last  member  of  the  chain  of  codistribu¬ 
tions,  J(N),  is  called  the  bottom  derived  system.  Since  =  /(N+1),  we  have  that 
dw  =  0  mod  J(N)  for  all  a?  €  J(N\  and  by  the  Frobenius  theorem,  the  bottom  de¬ 
rived  system  is  integrable.  That  is,  there  there  exist  functions  hi,...  ,hq  such  that 
J(N)  =  {dhi,. . .  ,dhq).  Solution  trajectories  of  I  are  then  constrained  to  lie  on  level 
surfaces  of  h. 

A  Pfaffian  system  is  said  to  be  nonholonomic  if  7 ^  is  a  proper  subset  of  7.  We 
will  only  work  with  systems  which  are  maximally  nonholonomic ,  that  is,  with  bottom 
derived  system  I ^  =  {0}. 

For  reference,  we  state  here  a  variant  of  the  familiar  Frobenius  theorem: 
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Theorem  1.  Frobenius  [1], 

Let  {w1,...  ,wp}  be  a  set  of  linearly  independent  one-forms,  and  fi,...  ,  /,  a  set  of 
functions  whose  differentials  are  linearly  independent  of  each  other  and  of  the  u>*  ’s. 

dw*  A  w1  A  •  •  •  A  u F  A  df  j  A  •  •  •  A  dfq  =  0 

for  i  =  1, . . .  ,p,  then  there  exist  coordinate  functions  Z\,. . .  ,  zv  and  coefficient  func¬ 
tions  dij,  b{j  such  that  the  one-forms  u>1  can  be  written  as: 

=  2  a'idzi  + 1]  M/i 

i=i  j=i 

The  proof  follows  from  the  standard  Frobenius  theorem  and  the  fact  that  the 
codistribution  {w1, . . .  ,u>p,  dfu . . . ,  d/,}  is  integrable. 

3.  Goursat  Normal  Forms 

A  great  deal  of  work  has  been  done  on  transforming  Pfaffian  systems  into  normal 
forms.  Of  special  interest  to  us  in  this  paper  is  the  Goursat  normal  form,  originally 
proposed  by  Goursat  for  Pfaffian  systems  of  codimension  two,  which  has  the  property 
that  all  its  integral  curves  can  be  expressed  in  terms  of  two  arbitrary  functions.  We 
will  also  examine  an  extended  Goursat  normal  form  for  Pfaffian  systems  of  codimen¬ 
sion  k  >  2,  for  which  solution  trajectories  can  be  expressed  in  terms  of  k  arbitrary 
functions. 

Definition  4.  Goursat  Normal  Form. 

A  codimension  two  Pfaffian  system  /  on  Rn  with  generators  of  the  form 
I  —  span{dzn  —  zn_idzi, . . .  ,dz3  —  z2dzi}. 
is  said  to  be  in  Goursat  Normal  Form. 

If  we  define  w1  (z)  =  dzn  —  zn-\dz\,...  ,wn_2(z)  =  dz$  —  z2dz\ ,  then  the  derived 
flag  of  I  is  given  by 

J<°>=  {uj1,  o>2,  ,  W"-3,  a,"-*} 

/<*>=  {u;1,  w2,  wn-3} 

/("-<)=  u>2} 

/("-3)  =  {u;1} 

/("-2)  =  {0} 

From  the  form  of  the  Pfaffian  system  in  the  z  coordinates,  it  follows  that  inte¬ 
gral  curves  of  the  system  are  unconstrained  in  their  Z\,zn  coordinates  alone.  Once 
Zi(t),zn(t)  are  specified  as  functions  of  some  parameter  t ,  the  other  coordinates  are 
determined  as  functions  of  zi(t),zn(t)  and  finitely  many  of  their  derivatives.  The 


GOURSAT  NORMAL  FORMS 


7 


following  classical  theorem  gives  necessary  and  sufficient  conditions  for  converting  a 
Pfaffian  system  into  Goursat  normal  form: 

Theorem  2.  Goursat  Normal  Form  [1]. 

A  Pfaffian  system  I  of  codimension  two  on  Rn  has  a  set  of  generators  which  are  in 
Goursat  normal  form  if  and  only  if  there  exists  a  basis  set  of  forms  {a1,...  ,an-2} 
for  I  and  a  one-form  it  satisfying  the  congruences: 

da'  =  — a,+1A?r  mod  a1,...,  a*  *=l,...,n  —  3 

dan~2  0  mod  I 

In  [10]  we  showed  that  the  Pfaffian  system  associated  with  the  system  of  a  car 
towing  n  trailers,  generated  b^^fBeTOnstfaints  that  each  axle  of  wheels  roll  without 
slipping,  satisfied  the  conditions  for  conversion  to  Goursat  normal  form.  This  was 
a  system  of  codimension  two,  corresponding  to  the  fact  that  the  linear  and  angular 
velocities  of  the  front  car  are  the  inputs,  and  thus  freely  specifiable. 

In  order  to  consider  Pfaffian  systems  associated  with  mobile  robots  such  as  the 
firetruck  [3,  9]  or  the  multi-steering  multi-trailer  system  [11]  we  need  to  work  with 

systems  of  codimension  greater  than  two.  We  have  the  following  definition: 

•  ' 

Definition  5.  Extended.  Goursat  Normal  Form. 

A  Pfaffian  system  I  on  Rn+m+1  of  codimension  m  -f  1  is  in  extended  Goursat  normal 
form  if  it  is  generated  by  n  constraints  of  the  form: 

I  =  span {dz\  —  Zi+idz°  :  j  =  1, . . .  , m; *  =  1, . . .  , Sj},  (2) 

We  note  that  this  is  a  direct  extension  of  the  Goursat  normal  form,  and  all  solution 
trajectories  of  (2)  are  determined  by  the  m  +  1  functions  z°(t),z\(t),. . .  , z™(t)  and 
their  derivatives  with  respect  to  the  parameter  t. 

There  are  conditions  due  to  Murray  [6]  for  converting  a  Pfaffian  system  to  extended 
Goursat  normal  form.  We  restate  and  prove  this  with  the  additional  condition  (cor¬ 
rection)  that  it  needs  to  be  integrable: 

Theorem  3.  Extended  Goursat  Normal  Form  [6]. 

Let  I  be  a  Pfaffian  system  of  codimension  m  +  1.  If  there  exists  a  set  of  generators 
{of  :  j  =  1,...  ,m;i  =  1,...  ,s,}  fori  and  an  integrable  one-form  it  such  that  for 
all  j, 

da1-  =  — A  it  mod  /(*>“*)  i  = 
da{.  ft  0  mod  I 

then  there  exists  a  set  of  coordinates  z  such  that  I  is  in 

I  =  { dz\  —  z{+ldz°  :  j  =  1, . . .  ,  m;  t  = 


!)•••  1 


(3) 


Goursat  normal  form, 

1 , .  . .  ,  Sj  } . 
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Proof.  If  the  Pfaffian  system  is  already  in  extended  Goursat  normal  form,  the  con¬ 
gruences  are  satisfied  with  n  =  dz°  and  the  basis  of  constraints  aj  =  dz\  -  z3i+1dz°. 

Now  assume  that  we  have  found  a  basis  of  constraints  for  I  which  satisfies  the 
congruences  (3).  It  is  easily  checked  that  this  basis  is  adapted  to  the  derived  flag, 
that  is: 

=  {a{  :j  =  1,...  ,m;  i  =1,...  ,s,-  —  jfc} 

We  will  now  construct  the  coordinates  z  which  comprise  the  Goursat  normal  form. 

Since  it  is  integrable,  any  first  integral  of  tt  can  be  used  for  the  coordinate  z°.  If 
necessary,  we  can  rescale  the  constraints  aj  so  that  the  congruences  (3)  are  satisfied 
with  dz° : 

d<4  =  -of+1  A  dz°  mod  /(*>“ *)  *  =  1, . . . ,  sj  -  1 
and  we  can  renumber  the  constraints  so  that  $i  >  s2  >  •  •  •  >  sm. 

Now  consider  the  last  nontrivial  derivedsystem,  The  one-forms  {a}, . . .  ,  aj  } 

form  a  basis  for  this  codistribution,  where  sj  =  s2  =  •••  =  sr,.  Prom  the  fact  that  * 

da{  =  A  dz°  mod  7^1_1), 

it  follows  that  the  one- forms  c*] , . . .  ,  oj1  satisfy  the  Frobenius  condition: 

da{  A  a\  A  •  •  •  A  a?  A  dz°  =  0 

and  thus,  by  the  Frobenius  theorem,  we  can  find  coordinates  z\ , . . .  ,  z\l  such  that 

a\  1  f  dz\  I 

;  =A  :  +  Bdz° 

a?  J  [  dzl1 

We  note  that  the  matrix  A  must  be  nonsingular,  since  the  oj’s  are  a  basis  for  7(4»_1) 
and  they  are  independent  of  dz°.  Therefore,  we  can  define  a  new  basis  dj  as: 

'  *i  1  r  i  r  l 

:  :=  A'1  :  =  :  +  (A~1B)dz° 


and  we  define  the  coordinates  z{  :=  (A^B)^  so  that  the  one-forms  a{  have  the  form 

dj  =  dz{  —  z%dz° 

for  j  =  1,...  ,ri. 

By  the  proof  of  the  standard  Goursat  theorem,  all  of  the  coordinates  in  the  jth 
tower  can  be  found  from  z{  and  z°,  thus  by  the  above  procedure,  we  have  effectively 
found  all  the  coordinates  in  the  first  ri  towers. 

To  find  the  coordinates  for  the  other  towers,  we  need  to  look  at  the  lowest  derived 
systems  in  which  they  appear.  The  coordinates  for  the  longest  towers  were  found 
first,  next  we  will  find  those  for  the  next-longest  tower(s). 
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Consider  the  smallest  integer  k  such  that  dim  >  kr^  more  towers  will 

appear  at  this  level.  We  know  that  a  basis  for  is 

f-1  zrl  pJT\  zrT\  1*1  + 1  /vri'’  f2\ 

where  of  =  dz3-  -  zf+1<fz°  for  j  =  1,...  ,rj,  as  found  in  the  first  step,  and  of  for 
j  =  t*i  +  1,. ..  ,r2  are  the  one-forms  we  started  with,  which  satisfy  the  congruences 
(3)  and  are  adapted  to  the  derived  flag.  The  lengths  of  these  towers  are  sr,+i  = 
•  •  •  sr,+rj  =  —  k  + 1.  For  notational  convenience,  we  will  define  zfa  :=  (zf , ...  ,  zf) 

for  j  —  1 , . . .  ,  t*i  .  .  . 

By  the  Goursat  congruences,  we  know  that  da{  =  —of  A  dz°  mod  I^Sl~  ’  for 
j  =  r*i  +  1, . ..  ,ri  +  r2,  thus  the  Frobenius  condition 

da{  A  of1  +1  A  •  •  •  A  a?"3  A  dz\  A  •  •  •  A  dzf  A  •  •  •  A  dzf1  A  •  •  •  A  dzf1  A  dz°  =  0 

is  satisfied  for  j  =  n  +  1, . . .  ,  rx  +  r2.  Using  the  Frobenius  theorem,  we  can  find  new 
coordinates,  z?",. . .  ,  z?"2 ,  such  that 

r  o?+i  i  r  dz?"  I  r  ^  1 


of1+ra 


dz?"2 


+  Bdz°  +  C 


Since  the  congruences  are  only  defined  up  to  mod  we  can  eliminate  the  last 

group  of  terms  (those  multiplied  by  the  matrix  C)  by  adding  in  the  appropriate 
multiples  of  of  =  dz{  -  zf+J dz°  for  j  =  1, ...  ,n  and  t  =  1, ..., fc.  This  will  change 
the  B  matrix,  and  we  will  be  left  with 

r  o?+l  i  r  dz?"  ] 


a?"2 


+  Bdz° 


dz?"2 


Again,  we  note  that  A  must  be  nonsingular  because  the  of ’s  are  linearly  independent 
mod  and  also  independent  of  dz°,  and  so  we  can  define 

’  a?"  1  f  dz?"  ' 

:=  A"  :  +  (A"B)dz°  =  :  +(A"B)dz° 

a?"3  J  dz?"3 

We  then  define  the  coordinates  zf  :=  (A"B)j  for  j  =  ri  +  1,...  ,ra  +  r2  so  that 
of  =  dz{  —  z$dz°.  Again,  by  the  standard  Goursat  theorem,  all  of  the  coordinates  in 
the  towers  ri  +T, . . .  ,rj  +  f2  are  now  defined. 

The  coordinates  for  the  rest  of  the  towers  are  defined  in  a  manner  exactly  analogous 
to  those  of  the  second-longest  tower. 

If  x  is  not  integrable,  then  we  cannot  use  the  Frobenius  theorem  to  find  the  coor¬ 
dinates.  In  the  special  case  where  si  >  s2,  that  is,  there  is  only  one  tower  which  is 
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longest,  it  can  be  shown  that  if  there  exists  any  n  which  satisfies  the  congruences, 
then  there  also  exists  an  integrable  t*  which  also  satisfies  the  congruences  (with  a 
rescaling  of  the  basis  forms),  see  [2,  6].  However,  if  si  =  S2,  or  there  are  at  least  two 
towers  which  are  longest,  this  is  no  longer  true.  Thus,  we  need  the  assumption  that 
7r  is  integrable.  □ 

If  I  can  be  converted  to  extended  Goursat  normal  form,  then  the  derived  flag  of  I 


has  the  structure: 

II  II 

• 

...  ... 

>  1  uli> 

•••»  <*?, 

• 

{*!, 

• 

•••» 

•  . 

/(*  1-2)  _ 
/(*  1-1)  = 
/(» l)  — 

W, 

{«!} 

{0} 

^2} 

where  the  forms  in  each  level  have  been  arranged  to  show  the  different  “towers”  which 
result.  The  superscripts  j  indicate  the  tower  to  which  each  form  belongs,  and  the 
subscripts  i  index  the  position  of  the  form  within  the  jth  tower.  There  are  Sj  forms 
in  the  jtk  tower.  An  algorithm  for  converting  systems  to  extended  Goursat  normal 
form  is  given  in  [2]. 

We  also  give  another  version  of  this  theorem,  which  is  easier  to  check,  since  it 
does  not  require  finding  a  basis  which  satisfies  the  congruences  but  only  one  which 
is  adapted  to  the  derived  flag.  We  also  give  the  proof  for  this  theorem,  since  only  a 
special  case  is  proved  in  [7]. 

Theorem  4.  Extended  Goursat  Normal  Form  [7]. 

A  Pfaffian  system  I  of  codimension  m  +  1  on  Rn+m+1  can  be  converted  to  Goursat 
normal  form  if  and  only  if  there  exists  a  one-form  v  such  that  {I^k\  7r}  is  integrable 
for  k  =  0, . . .  ,  N  —  1 . 

Proof  The  only  if  part  is  easily  shown  by  taking  jt  =  dz°  and  noting  that 

/(*)  =  {dz{  —  z£_jdz°  :  j  —  1,...  ,m;t  =  +  1,...  ,sj} 

{/W,tt}  =  {dz{,dz°  :  j  =  1,...  ,m;i=  + 1,...  ,Sj}, 

which  is  integrable. 

For  the  if  part,  assume  that  such  a  ir  exists.  First,  we  find  the  derived  flag  of  the 
system,  I  D  D  •••  D  /M  =  {0}.  We  will  iteratively  construct  a  basis 

which  is  adapted  to  the  derived  flag  and  which  satisfies  the  Goursat  congruences  (3). 

We  claim  that  the  lengths  of  each  tower  are  determined  from  the  dimensions  of 
the  derived  flag.  Indeed,  the  longest  tower  of  forms  has  length  Sj.  If  the  dimension 


GOURSAT  NORMAL  FORMS 


11 


of  is  r1}  then  there  are  rj  towers  which  each  have  length  si\  and  we  have 

sj  =  s2  =  •  ••  =  sri.  Now,  if  the  dimension  of  is  2ri  +  r2,  then  there  are  r2 

towers  with  length  si  —  1,  and  we  have  sr,+j  =  •  •  •  =  3ri+ra  =  «i  —  1.  We  find  each 
Sj  similarly. 

We  note  here  that  a  jt  which  satisfies  the  conditions  must  be  in  the  complement 
of  /,  for  if  it  were  in  /,  then  {/,*•}  integrable  means  that  I  is  integrable,  and  this 
contradicts  our  assumption  that  I  is  maximally  nonholonomic,  that  is  I ^  =  {0}  for 
some  N. 

Consider  the  last  nontrivial  derived  system,  Let  {a}, . . .  , Qj1 }  be  a  basis 

for  The  definition  of  the  derived  flag,  specifically  1^  =  {0},  implies  that 

da{  0  mod  j  =  l,...,ri  (4) 

Also,  the  assumption  that  {I^k\ir}  is  integrable  gives  us 

da{  =  0  mod  jr}  j  —  (5) 

combining  equations  (4)  and  (5),  we  have  that 

da{  =  7r  A  mod  j  —  (6) 

for  some  ft  ^  0  mod 

Now,  we  also  have  from  the  definition  of  the  derived  flag  that 

da{  =  0  mod  j  =  1, . . .  ,ri 

which  combined  with  (6)  gives  us  that  ft  is  in  f(ei~2K 
Claim,  ft,... ,  ft1  are  linearly  independent  mod 

Proof  of  Claim.  The  proof  is  by  contradiction.  Suppose  there  exists  some  combination 
of  the  ft's,  say 

0  =  bift  H - K  6r,  ft1  =0  mod 

with  not  all  of  the  bfs  equal  to  zero.  Consider  a  =  bia\  +  •••  +  bTl  ap .  We  must 
have  a  ^  0  because  the  a{  are  a  basis  for  The  exterior  derivative  of  a  can  be 

found  by  the  product  rule, 

da  =  bjda{  +  53  dfy  A  aj 
i=i  i=i 

=  j£,bj(KAft)  mod  I*'1-1* 

i=i 

=  T  A  (T.bjft)  mod 

i= i 

=  0  mod 
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which  implies  that  a  is  in  I^K  However,  this  contradicts  our  assumption  that 
J(*i)  =  {0}.  Thus  the  ft' s  are  linearly  independent  mod  and  the  claim  is 

proven. 

Define  Oj  :=  ft  for  j  =  1, . . .  ,  Note  that  these  basis  elements  satisfy  the  first 
level  of  Goursat  congruences,  that  is: 

da{  =  A  7r  mod  j  =  1, . . .  ,n 


If  the  dimension  of  is  greater  than  2rj,  then  choose  one-forms  oJ1+1, 

such  that  that 


J/v1  r/1  n/1  nTi  /vri+1  rJ l+rs\ 

is  a  basis  for  /(4»-2). 

For  the  induction  step,  we  assume  that  we  have  a  basis  for  /W, 
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aj1+rj 


{®1  »  *  *  •  »  5  »  •  ♦  •  »  )  ®|fcc  } 

which  satisfies  the  Goursat  congruences  up  to  this  level: 

dot{  =  — a*+i  A  7r  mod  /(*■»“*)  j  =  l,...,c;  £  =  !,.. 


Note  we  have  assumed  that  c  towers  of  forms  have  appeared  in  /W.  Consider  only 
the  last  form  in  each  tower  that  appears  in  7(,),  that  is  a*k.,j  =  1, . . .  ,c.  By  the 

construction  of  this  basis  (or  from  the  Goursat  congruences),  we  have  that  cdk  is  in 
7^  but  is  not  in  I^t+1\  thus 


da{.  ^  0  mod  7^  j  =  1, . . .  ,  c 


The  assumption  that  {/W,7r}  is  integrable  assures  us  that 

da{.  =  0  mod  {/<*), r}  j  =  1,...  ,c 
thus  we  have  that  da{.  must  be  a  multiple  of  n  mod 

da{.=irAft  mod  /W  j  — 

for  some  ft  =£  0  mod  7W.  We  also  have,  from  the  fact  that  akk  is  in  1^  and  the 
definition  of  the  derived  flag,  that  1 

da?k.  s  0  mod  7(,-1)  j  =  1, . . .  ,  c 

which  implies  that  ft  €  7^_1h  By  a  similar  argument  to  the  claim  above,  we  can 
show  that  the  ft's  are  independent  mod  7^.  We  define  aj>+1  =  ft,  and  thus 

{ai » •  •  •  >aJtj+nan***  ?a*2+i»***  »<*!»•••  >alc+i} 
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forms  part  of  a  basis  of  7b-1).  If  the  dimension  of  is  greater  than  &i  +  &2  H - I- 

ke  +  c,  then  we  complete  the  basis  of  with  any  linearly  independent  one-forms 
aj+1 , . . .  ,  aj+rc  such  that 


o 


iv 


/v2 

'aW 


.a 


IV 


rtc  rt0*1 


is  a  basis  for  7b-1). 

Repeated  application  of  this  procedure  will  construct  a  basis  for  7  which  is  not 
only  adapted  to  the  derived  flag,  but  also  satisfies  the  Goursat  congruences. 

We  note  that  by  assumption,  ir  is  integrable  mod  the  last  nontrivial  derived  system, 
Looking  at  the  congruences  (3),  we  see  that  any  integrable  one-form  n'  which 
is  congruent  to  ir  up  to  a  scaling  factor, 

7 r'  =  dt  =  fit  mod  7(*1-1) 

will  satisfy  the  same  set  of  congruences  up  to  a  rescaling  of  the  constraint  basis  by 
multiples  of  this  factor  /.  0 


4.  Pfaffian  Systems  generated  by  Mobile  Robots 

We  now  consider  some  multi-steering  mobile  robot  systems  and  show  that  the 
Pfaffian  systems  generated  by  the  constraints  that  the  wheels  roll  without  slipping 
satisfy  the  extended  Goursat  conditions. 

Example  1.  Firetruck. 

Consider  the  example  of  a  firetruck  [3].  There  are  two  steering  wheels  in  the  system: 
one  at  the  front,  for  the  driver,  and  another  at  the  rear,  for  the  tiller.  We  model 


FIGURE  1 .  A  sketch  of  the  Firetruck,  with  steering  wheels  on  the  front 
and  back  axles. 

the  firetruck  as  a  Pfaffian  system  generated  by  the  constraints  that  the  three  pairs  of 
wheels  roll  without  slipping: 

7  =  K,a\a2} 
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In  coordinates,  the  constraints  can  be  written  as: 

a0  =  sin(0i  +  <f>i)dxo  —  cos(0i  +  <l>i)dyo 
a1  =  sinOidxi—cosOidyi 

a2  =  sin($2  +  ^2)^2  —  cos(#2  +  <f>2)dy2 

Since  the  coordinates  xo,yo,x2,y2  are  determined  by  xi,yi,0t,#2,  we  can  parame¬ 
terize  the  state  space  by  q  =  {xi,yi,0i,02,<l>i,<f>2}.  A  complement  to  the  system  is 
{dx!,d<f>i,d<j>2}.  It  can  be  shown  that  the  constraints  are  adapted  to  the  derived  flag, 

I-  {a0,©1, a2} 

/<»>=  {a1} 

J<2>  =  {0} 

and  that  the  systems  {/,dij}  and  {I^\dx\}  are  integrable.  By  Theorem  4,  the 
system  can  be  converted  into  extended  Goursat  normal  form.  We  refer  the  reader 
to  [3,  9]  for  the  coordinate  transformation  (into  chained  form,  which  is  the  dual  of 
Goursat  form),  and  methods  for  steering  the  system. 

A  multi-steering  trailer  system  was  examined  in  (11],  and  it  was  shown  how  to 
transform  such  a  system  into  chained  form  (which  is  the  dual  of  Goursat  normal 
form)  using  dynamic  state  feedback.  That  is,  states  were  added  to  the  system  and 
this  augmented  system  was  transformed  into  the  dual  of  Goursat  normal  form.  In  this 
section,  we  show  that  this  augmentation  is  not  always  necessary;  some  configurations 
of  the  multi-steering  trailer  system  can  be  transformed  into  Goursat  form  using  only 
static  state  feedback. 

We  concentrate  on  the  specific  case  of  a  5-axle  system  with  two  steering  wheels. 
Assuming  that  the  first  axle  is  steerable,  there  are  four  possible  positions  for  the 
second  steering  wheel.  We  examine  these  four  cases  to  show  the  variety  of  results 
that  can  be  obtained  from  a  relatively  simple  system  configuration. 

Example  2.  5-axle ,  1-2  steering  trailer  system. 

Consider  the  case  of  a  5-axle  trailer  system  with  the  first  two  axles  steerable.  A 
sketch  of  this  system  is  shown  in  Figure  2.  The  state  space  is  parameterized  by 
q '  =  {x5,  t/5, 05, 64, 03, 02,  the  x,y  position  of  the  last  axle,  the  angle  of  each 

hitch  0„  and  the  extra  steering  angle  <f>.  We  note  that  the  z,-,  y,  positions  of  the  other 
axles  can  be  written  in  terms  of  the  state  variables  q. 

The  constraints  are  that  each  axle  roll  without  slipping: 

a'  =  sin  Oidxi  —  cos  fljdy,  t  =  1,3, 4, 5 
a2  =  sin  <t>dx2  —  cos  <f>dy2 

The  Pfaffian  system  is  I  =  {a1, a2, a3, a4, a5},  a  complement  to  this  system  is  given 
by  {d<j>,d9i,dxs},  and  /  has  codimension  three. 
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FIGURE  2.  A  5-axle  trailer  system  with  the  first  two  axles  steerable. 

This  basis  is  adapted  to  the  derived  flag, 

7  =  {a1,  a2,  a3,  a4,  a5} 
jW  =  {a3,  a4,  a5} 

J(2)  =  {a4,  a5} 

J(3)  =  {a5} 

J(4>  =  {0} 

and  the  congruences  are. satisfied  with  this  basis  as  well: 

da 5  =  Cs(q)  a4  A  dx $  mod  7^ 

da 4  =  c4{q)  a 3  A  dx6  mod  7(2) 

da 3  =  c3(g)  a2  A  dx5  mod  7(1) 

da2  =  Ci(q)  d(f>  A  dx5  mod  I 

da1  =  Ci(q)dOiAdxs  mod  7 

We  note  here  that  by  a  simple  rescaling  of  the  basis,  the  functions  Ci(q)  can  be 
eliminated  to  get  the  Goursat  congruences  (3)  exactly.  This  is  done  as  follows.  First 
define  a4  :=  -c5(q)a4  to  get  da5  s  -a4  Adxs  mod  7<3>.  Then,  taking  the  derivative 

of  a4,  we  see  that 

da4  =  - cs(q)da 4  -  dc6(q)  A  a4 

=  —Cs(q)c^(q)a3  A  dxs  mod  7^ 

since  a 4  is  in  7(2).  Defining  a3  :=  cs{q)c4(q)a 3,  we  have  da4  =  -a3  A  dxs  mod  7(2). 
The  other  constraints  are  scaled  similarly.  For  the  rest  of  this  paper,  we  will  assert 
that  the  Goursat  congruences  are  satisfied  if  we  have  the  modified  congruences 

daj  =  cj+i(9)  <4+i  A  T  mod  *  =  !>  •  •  •  > -  1 

dcd  #  0  mod  7 

instead  of  the  original  congruences  (3). 
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Example  3.  5-axle ,  1-5  steering  trailer  system. 

Now  consider  the  same  5-axle  system  as  in  the  previous  example  with  the  first  and 
fifth  (last)  axles  steerable,  as  sketched  in  Figure  3.  The  configuration  space  can  be 


Figure  3.  A  5-axle  trailer  system  with  steering  wheels  on  the  first 
and  last  axles. 


parameterized  by:  q  =  {£4,  t/4,  #5, #3;  the  x,y  position  of  the  second-to- 

last  axle,  the  angle  of  each  hitch,  and  the  angles  of  the  steering  wheels.  As  before, 
the  Xi ,  y,  positions  of  the  other  axles  can  be  written  in  terms  of  the  coordinates  q. 
The  constraints  are  that  each  axle  roll  without  slipping, 

a’  =  sin  Oidii  —  cos  0,dy,  *  =  1 , 2, 3, 4 

a5  =  sin  <j>dx5  —  cos  <j>dy5 


The  Pfaffian  system  is  7  =  {a1,©2, a3, a4, a5}.  A  complement  to  this  system  is  given 
by  {d<j>,dOi,dx4}-,  I  has  codimension  three.  This  basis  is  adapted  to  the  derived  flag, 


1  = 
/<!>  = 
/(2)  = 
/(3)  _ 

7<4>  = 


{a1,  a2,  a3,  a4,  a5} 
{a2,  a3,  a4} 

{a3,  a4} 

{a4} 

{0} 


The  congruences  are  satisfied  with  this  basis  as  well, 


da *  ss  c<(g)  q3  A  dx4  mod  1^ 
da 3  =  03(g)  o2  A  dx4  mod  I ^ 
da 2  =  02(g)  a1  A  dx4  mod 
da 1  =  Cj(g)  dOi  A  dx  mod  7 
da 5  =  c5(g)  d<f>  A  dx4  mod  7 
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and  the  Pfaffian  system  can  be  converted  to  extended  Goursat  normal  form. 

Example  4.  5-axle ,  1-4  steering  trailer  system. 

We  now  consider  the  5-axle  system  with  the  first  and  fourth  axles  steerable,  as 
sketched  in  Figure  4.  The  configuration  space  can  be  parameterized  by  the  x,y 


FIGURE  4.  A  5-axle  trailer  system  with  the  first  and  fourth  axles  steerable. 

position  of  the  third  axle,  the  hitch  angles  0,,  and  the  steering  angle  of  the  third 
axle  <f>.  We  let  q  represent  the  state,  q  =  {13,  jte, $s,  #4,  03, 02,  0i,  4>}-  The  other  im¬ 
positions  can  be  written  in  terms  of  q. 

The  constraints  are  that  each  axle  rolls  without  slipping: 

a’  =  sin  Oidxi  —  cos  i  =  1, 2, 3, 5 

a4  =  sin  <f>dx4  —  cos  4>dy4 

The  Pfaffian  system  is  thus  7  =  {a1,  a2,  a3,  a4,  a5}  and  a  complement  to  this  system 
is:  {d<f>, d0i,<fi3}.  This  basis  is  adapted  to  the  the  derived  flag, 

I  =  {a1,  a2,  a3,  a4,  a5} 

7(1)  =  {a2,  a3,  a5} 

=  {a3} 

7<3>  =  {0} 

and  also  satisfies  the  Goursat  congruences: 

da3  =  Cz(q )  o?  A  dx 3  mod 

da 2  =  C2(q)  a1  A  dx$  mod  7^ 

da1  =  Ci(q)  d6i  A  dx3  mod  7 
das  =  c5(<7)  a*  A  dx 3  mod  7^ 

da4  =  04(9)  d<f>  A  dx 3  mod  7 

and  thus  the  Pfaffian  system  can  be  converted  to  extended  Goursat  normal  form. 

Example  5.  5-axle,  1-3  steering. 

The  final  instance  of  the  5-axle  trailer  system  has  the  first  and  third  axles  steerable,  as 
sketched  in  Figure  5.  The  state  space  is  parameterized  by :  q=  {is,  1/5,  #5, 04, 03, 02, 0i,  <£}, 
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FIGURE  5.  A  5-axle  trailer  system  with  the  first  and  third  axles  steer¬ 
able.  This  is  the  only  configuration  of  the  5-axle  system  with  two 
steering  wheels  which  does  not  satisfy  the  conditions  for  converting  to 
extended  Goursat  normal  form. 

and  the  other  sc,-,  y,  can  be  written  in  terms  of  q. 

The  constraints  are  that  each  axle  roll  without  slipping: 

q‘  =  sin  6{dxi  —  cos  0,dy,  i  =  1,2, 4, 5 
a3  =  sin  (j>dx 3  —  cos  <f>dy3 

The  Pfaffian  system  is  I  =  {a1,  a2, a3,  a4, a8},  and  a  complement  to  the  system  is 
given  by  {d<f>,d0udx3}. 

This  basis  is  adapted  to  the  derived  flag, 

I  =  {a1,  a2,  o;3,  a4,  a5} 

/W  =  {a2,  a4,  a5} 

/(2)  =  {a5} 

7(3)  =  {0} 

however,  the  congruences  are  not  satisfied: 

da 5  =  c5(q)  a4  A  dx3 
da4  =  c4(q)  a3  A  dx3 
da2  =  C2(q)  a1  A  dx3  -f  £2(9)  a1  A  a3 
da3  =  c3(q)  d<{>  A  dx 3 
da1  =  Cx(q)  d0}  A  dx3 

lit  order  to  have  {J^2\ir}  integrable,  we  must  choose  v  =  dx3  (mod  {o4, o5}).  This 
will  also  give  us  {7°,  x}  integrable,  but  {71,  «■}  is  not  integrable.  Thus,  this  system 
does  not  satisfy  the  conditions  for  conversion  to  extended  Goursat  normal  form.  We 
will  return  to  this  example  in  the  next  section. 

5.  Prolongations 

If  a  Pfaffian  system  7  of  codimension  k  satisfies  the  necessary  and  sufficient  con¬ 
ditions  for  converting  into  extended  Goursat  form,  then  its  solution  trajectories  are 


mod  7<2> 
mod  7^ 
mod  7^ 
mod  7 
mod  7 
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determined  by  k  arbitrary  functions.  However,  even  if  a  system  cannot  be  trans¬ 
formed  into  Goursat  form,  its  solution  trajectories  may  still  have  this  property.  If  so, 
then  we  say  that  7  is  absolutely  equivalent  in  the  sense  of  Cartan  to  the  trivial  system 
(the  system  with  no  constraints)  on  R*.  Although  we  will  not  examine  the  concept  of 
absolute  equivalence  in  its  full  generality,  we  will  give  some  sufficient  conditions  for  a 
Pfaifian  system  to  have  a  prolongation  which  can  be  converted  to  Goursat  form,  and 
thus  the  solution  trajectories  of  7  are  determined  by  k  independent  functions. 

A  general  type  of  prolongation  which  preserves  a  one  to  one  correspondence  be¬ 
tween  solution  trajectories  of  the  original  and  prolonged  system  is  a  Cartan  prolon¬ 
gation. 

Definition  6.  Cartan  Prolongation. 

Let  7  be  a  Pfaffian  system  on  a  manifold  M.  A  system  J  on  M  xR*  js  a  Cartan 
prolongation  of  7  if  : 

(1) ir*(/)CJ 

(2)  For  every  solution  curve  c  :  (— e,  c)  — ►  Af  of  7  there  exists  a  unique  solution 
curve  c  :  (-e, e)  -+  M  x  Rp  of  J  with  no  c-c. 

If  7  is  equipped  with  a  given  independence  condition  r,  then  we  also  require  that  n'r 
be  the  independence  condition  for  7. 

A  canonical  way  to  prolong  a  system  with  independence  condition  dt  is  to  take  an 
integrable  one-form  dij  in  the  complement  of  7,  and  augment  7  with  the  additional 
form  drf  -  ydt,  where  y  is  a  new  coordinate  on  R.  In  effect,  this  adds  the  derivative 
of  rj  (with  respect  to  the  independence  condition)  as  a  state  variable.  As  long  as  all 
solution  trajectories  are  “smooth  enough”  (we  will  assume  C°°),  there  will  be  a  one- 
one  correspondence  between  solution  trajectories  of  the  original  and  the  prolonged 
system. 

We  consider  a  special  type  of  Cartan  prolongation  which  consists  of  many  of  these 
canonical  prolongations. 

Definition  7.  Prolongation  by  differentiation. 

Let  7  be  a  Pfaffian  system  of  codimension  m  + 1  on  Rn+m+1  with  coordinates  (z,  v,t) 
for  which  dt  is  an  independence  condition  and  {dui, . . .  ,  dvm,  dt }  forms  a  complement. 
Let  ,bm  be  nonnegative  integers  and  let  b  denote  their  sum.  The  system  7 

augmented  by 

dv i  —  v\dt,  . . .  ,  dvj1-1  —  v\*dt, 

dv 2  —  v\dt ,  . . .  ,  dv*2~l  —  v^dt 

•  • 

•  • 

•  • 

dvm  -  vlmdt,  ,  . . .  ,  dvfr-1  -  vb"dt, 

is  called  a  prolongation  by  differentiation  of  7.  The  augmented  system  is  defined  on 
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We  note  here  that  since  the  original  system  and  independence  condition  corre¬ 
sponded  to  a  set  of  first  order  ordinary  differential  equations,  the  prolonged  system 
has  the  same  independence  condition  and  also  corresponds  to  a  set  of  first  order 
ordinary  differential  equations. 

We  can  now  give  sufficient  conditions  for  a  Pfaffian  system  to  have  a  prolongation 
which  is  equivalent  to  extended  Goursat  form. 

Theorem  5.  Conversion  to  Goursat  form  using  prolongation  by  differen¬ 
tiation. 

Consider  a  Pfaffian  system,  I  =  {a1,...  ,aB}  on  RB+m+1  with  independence  con¬ 
dition  dz°  and  complement  {dt>i, . . .  ,dvm,dz0} .  If  there  exists  a  list  of  integers 
bi,...  , bm  such  that  the  prolonged  system 

/  =  {  a1,...  ,orB,dt>i  —  vjdz0,...  —  v\*-dz°y 

...  ,dvm-  v'mdz°, ...  ,  dvb™~'  -  v^fdz0} 

satisfies  the  condition  that  {Pk\dz0}  is  integrable  for  all  k,  then  I  can  be  transformed 
to  extended  Goursat  normal  form  using  a  prolongation  by  differentiation. 

Proof.  The  proof  is  by  application  of  Theorem  4  to  the  prolonged  system  I.  □ 

Although  this  is  a  very  specific  form  of  prolongation  of  a  Pfaffian  system,  and  the 
conditions  of  the  theorem  must  be  checked  in  a  specific  coordinate  system  with  a  given 
independence  condition,  there  do  exist  practical  systems  which  can  be  converted  into 
extended  Goursat  form  using  this  type  of  prolongation. 


Example  6.  S-axle,  1-3  steering ,  revisited. 

We  return  to  the  5-axle  trailer  system  with  the  first  and  third  axles  steerable,  which 
did  not  satisfy  the  conditions  for  conversion  to  extended  Goursat  form. 

Recall  that  the  derived  flag  was  of  the  form: 


JW  =  {a\  a4,  a5} 

7(2)  =  {a5} 

7<3>  =  {0} 


but  that  the  congruences  were  not  satisfied: 


da5  =  cs(g)  a4  A  dx 3  mod  7^ 

da4  =  04(9)  or3  A  dx 3  mod  7^ 

da 2  =  Cilq)  a1  A  dx3  +  ki(q)  a1  A  a3  mod  7^ 

da3  =  c${q )  d<f>  A  dx3  mod  7 

da 1  =  Ci(q)  d0\  A  dx3  mod  7 


We  can  look  at  the  equations  for  the  exterior  derivatives  of  the  constraints  to  see 
if  after  prolongation,  the  augmented  Pfaffian  system  will  satisfy  the  conditions  for 
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conversion  to  extended  Goursat  form.  We  see  that  it  =  dx$  will  give  us  it) 
integrable.  However,  {I^\dx 5}  is  not  integrable  since  da 2  has  a  term  k2(q)  a1  A  a3. 
If  either  a1  or  a3  could  be  added  to  1^\  that  term  would  no  longer  cause  a  problem. 
We  note  that  if  a 1  were  added  to  then  da2  =  0  mod  7^  and  we  will  still  have 
the  same  problem,  except  now  with  I^2K  If  we  can  somehow  add  a3  to  I^\  it  appears 
that  the  conditions  of  Theorem  4  will  be  satisfied  (the  only  thing  remaining  to  be 
checked  is  that  da3  =  0  mod  a2,a3,a4,a5,dx5.) 

We  prolong  7  by  differentiation,  and  augment  it  by  the  additional  form  v  =  d<f> — 
vdx3.  The  derived  flag  of  the  augmented  system  is: 

7=  {a1,  a2,  u>,  a3,  a4,  a5} 
fi1)  =  {or2,  a3,  a4,  or5} 

.  7<2>  =  {a4,  a5} 

7<3>  =  {o5} 

7<4>  =  {0} 

and  the  systems  {ftk\  dx3}  are  integrable  for  all  k,  as  can  be  seen  from  looking  at 
the  Goursat  congruences, 

do5  =  cs(q)a4Adi3  mod  7^ 
da4  =  c4(q)  o3  A  dx 3  mod  1^ 

da3  =  c3(q)  du>  A  <fx3  mod  7^^ 

du>  =  Cu(q)  dv  A  dxz  mod  I 

da 2  =  c2(q)  o1  A  dx3  mod  7^1 

da 1  =  ci  (q)  d0i  A  dx3  mod  I 

Thus,  the  prolonged  system  7  can  be  converted  into  extended  Goursat  normal  form. 
This  extension  is  shown  in  [11],  as  are  methods  for  steering  this  type  of  system. 

6.  Control  Systems 

As  we  mentioned  in  the  introduction,  control  systems  are  a  special  type  of  Pfaffian 
system,  and  therefore  all  of  the  results  presented  thus  far  can  be  specialized  to  control 
systems.  Most  of  the  previous  work  analyzing  nonlinear  control  systems  has  been 
from  the  vector  field  point  of  view,  taking  a  system  x  =  f(x)  +gi{x)v.i  -+-•••  gm(x)um, 
and  looking  at  properties  of  the  vector  fields  /,</,.  The  Pfaffian  systems  formulation 
is  the  dual  of  this. 

Definition  8.  Control  System. 

A  control  system  x  =  f(x,  u)  with  the  state  x  €  Rn,  the  input  u  €  Rm,  and  the 
derivative  of  the  state  taken  with  respect  to  time  t  €  R,  generates  a  Pfaffian  system 
7  on  Rn+m+1 

7  =  {dxi  —  /'(x,  u)dt :  i  =  1, . . .  ,  n} 


(8) 
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with  complement  {dui, . . .  ,dum,dt).  The  natural  independence  condition  to  choose 
is  dt,  since  we  want  dt  ^  0  along  all  solution  trajectories  of  the  system. 

Any  Pfaffian  system  I  of  codimension  m  + 1  on  Rn+m+1  with  coordinates  (x,u,t) 
can  be  called  a  control  system  if  it  has  a  set  of  generators  of  the  form  (8). 

Brunovsky  showed  that  any  controllable  linear  system  x  =  Ax  +  Bu  with  x  €  R", 
u  €  Rm  can  be  converted  to  a  “canonical”  form  given  by 

i\  =  Uj  x\  =  «2 

x\  =  x\  xjs=xl 


'■  xl2  =  xl2-i 
xh  =  4,-1 

with  n  =  ki  +  •••  +  km.  A  control  system  is  said  to  be  linearizable  if  and  only 
if  it  can  be  converted  to  Brunovsky  form  using  a  nonlinear  coordinate  change  and 
state  feedback.  Since  Brunovsky  linear  form  for  a  control  system  is  a  special  case 
of  extended  Goursat  normal  form  (2)  with  dz°  =  dt  and  z*1+1  =  tij,  the  theorems 
for  transforming  to  Goursat  form  can  be  specialized  to  give  conditions  for  exact 
linearization. 

Theorem  6.  Exact  Linearization  [5]. 

If  a  control  system  I  defined  on  Rn+m+1  has  a  set  of  generators  {of  :  j  =  1, . . .  ,  m;  *  = 
1, . . .  ,  sj }  such  that  for  all  j, 

da{  =  —  of+1  A  dt  mod  /(*>“*)  t  =  1, . . .  ,Sj  —  1  .  . 

dcd  £  0  mod  I 

Sj  7- 

then  there  exists  a  set  of  coordinates  z  such  that  I  is  in  Brunovsky  normal  form, 

I  =  {dt\  —  z3i+1dt  :  j  =  1,...  ,m;i  =  1,...  ,s,}. 

An  algorithm  for  converting  systems  to  Brunovsky  normal  form  is  also  given  in  [5], 
and  it  is  shown  that  if  the  control  system  is  time-invariant  and  affine  in  the  inputs, 
then  the  resulting  feedback  transformation  is  also  autonomous  and  input-affine. 

The  control  systems  version  of  Theorem  4  is  given  by 

Theorem  7.  Exact  Linearization  [8]. 

A  control  system  I  can  be  converted  to  linear  form  if  and  only  if  {I^k\dt}  is  integrable 
for  every  k. 

By  way  of  example,  we  examine  a  control  system  which  is  not  linearizable  but  can 
be  converted  to  Goursat  normal  form.  The  transformation  scales  time  by  a  function 
of  one  of  the  states. 


=  u 
X?  as  X 


=  XL-1 


(9) 
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Example  7.  Goursat  normal  form  for  a  control  system. 

Consider  the  single-input  control  system  [4], 

*i  =  *2  +  *32 

X2  =  #3 

i3  =  U 

This  control  system  generates  a  Pfaffian  system, 

I  =  {dx  i  —  (z2  +  xs2)dt,  dx  2  —  x$dt,  dx  3  —  udt}  (11) 

I  is  of  codimension  two  on  R5  with  coordinates  (x,u,<).  The  derived  flag  of  I  is 

I  ={a\a2,a3} 

JW  ={or1,a2} 

I™  ={a1 } 

/(3)  ={0} 

where  the  one-forms  adapted  to  the  derived  flag  are  given  by 

a1  =  dx  1  —  2  x3dx2  +  (x32  —  *2  )dt 
or2  =  dx  2  —  X3  dt 
a3  =  dx3  —  udt 

Note  that  this  is  not  the  basis  of  (11)  which  generated  7.  Since  {I^2\dt}  is  not 
integrable,  the  system  is  not  feedback  linearizable  by  Theorem  7. 

We  find  however  that  the  Goursat  congruences  (1)  are  satisfied,  for  n  —  dr  = 
dt  —  2  dx3: 

da 1  =  a2  A  dr 

da2  =  c(u)  a3  A  dr 

da3  =  c(tt)  du  A  dr  mod  a3 

Thus,  there  does  exist  a  transformation  $(x,u,<)  =  (z,v,r)  to  Goursat  normal  form, 
which  is  given  by 


T 

=  t  ~ 

-  2x3 

V 

u 

T 

-  2u 

=  x3 

*2 

=  X2 

-*32 

*3 

=  Xi 

.  2 
—  2x2X3  +  -X3' 
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and  it  is  easily  checked  that 

dz\ 

17  =  V 

dZ2 

17  =  *' 

dzz 

17  =  22 

This  formulation  of  the  system  makes  it  simple  to  analyze  the  possible  trajectories 
of  the  system  in  state/time  space.  Also,  a  controller  could  be  designed  using  linear 
techniques  which  would  be  valid  in  certain  regions  of  the  state  space. 

Converting  control  systems  to  Goursat  normal  form  may  not  be  the  most  useful 
thing  to  do.  However,  linearizing  control  system  using  dynamic  extension  is  a  problem 
that  has  been  studied  extensively.  A  dynamic  extension  of  a  control  system  is  an 
augmented  system  with  integrators  added  to  the  inputs;  for  example,  a  simple  first- 
order  dynamic  extension  is  given  by: 

x  =  /(x,u) 
iik  -  v 

where  an  integrator  is  added  to  the  ktk  input  channel. 

The  prolongation  by  differentiation  which  we  defined  in  Section  5  is  exactly  the 
dual  of  dynamic  extension  in  the  language  of  forms.  Thus,  we  have  the  control 
systems  version  of  Theorem  5: 

Theorem  8.  Linearization  by  dynamic  extension. 

Consider  a  control  system  I  on  Rn+m+1  with  coordinates  (x,u,t),  independence  con¬ 
dition  dt,  and  complement  {dui,...  ,dum,dt}.  If  there  exists  a  prolongation  by  dif¬ 
ferentiation  of  dimension  b  =  bi  +  •  •  *  +  bm  such  that  the  augmented  system 

/  =  {  a' =  dii  -  f'(x,u)dt :  i  =  l,...,n; 

0j  =  duj~l  ~  ujdt :  j  =  1,..  •  ,m;  k  =  0,. . .  ,6,} 

on  Rn+m+6+1  satisfies  the  condition  { /M,  dt}  is  integrable  for  every  k,  then  the  orig¬ 
inal  system  I  is  linearizable  by  dynamic  extension. 

Proof.  Apply  Theorem  7  to  the  extended  system  /.  □ 

This  theorem  is  similar  to  the  one  stated  by  Charlet,  Levine,  and  Marino  [4]  which 
also  gave  sufficient  conditions  for  linearizing  systems  by  dynamic  extension.  Their 
conditions  also  relied  on  the  existence  of  some  integers  6,  which  are  the  number  of 
integrators  added  to  the  ?k  input  channel.  However,  the  existence  of  a  dynamic 
extension  of  order  b  =  (6j,...  , bm)  which  is  linearizable  does  not  imply  that  the 
conditions  of  their  theorem  are  satisfied  for  that  6;  whereas  if  there  exists  a  dynamic 
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extension  of  order  b  =  (b\, ...  ,  bm)  which  can  be  linearized,  our  conditions  will  always 
be  satisfied  for  that  b. 

We  present  a  simple  example  to  show  how  the  theorem  can  be  applied  to  linearize 
control  systems  using  dynamic  extension. 

Example  8.  A  control  system  satisfying  the  conditions  of  Theorem  8  but 
not  the  conditions  of  [4]. 

Consider  a  4-state,  2-input  control  system: 

*1  =  X2  +  X3U2 

®2  =  X3  -f  Xi«2 

X3  =  tii  +  X2U2 
x4  =  u2 

The  corresponding  Pfaffian  system  on  R7  is 

1  =  {dx  1  —  ( x2  -I-  x3u2)dt,  dx2  —  (X3  +  xi u2)dt,  dx3  —  (ui  +  X2 u2)dt,  dx4  —  u2dt} 

(12) 

with  independence  condition  dt  and  complement  {dui,du2,dt}.  The  derived  flag  has 
the  form: 

I  =  {oS  Qf2)  a3,  a4} 

1^)  =  {a1,  a4} 

If2)  =  {0} 

The  one-forms  a'  which  are  adapted  to  the  derived  flag  are  not  the  same  as  those  of 
(12)  which  generated  I, 

a 1  =  dx  1  —  X2  dt 
a 2  =  <fx2  —  (X3  +  xiU2)dt 

a3  —  dx  3  —  (ui  +  X2U2)dt 
a4  =  dx4  +  x3di2  —  x3dt 
The  structure  equations  are  fairly  simple  to  find, 
doi1  —  —a2  A  dt 
da2  =  —Xidu2  A  dt 
da3  =  —  du\  A  dt  —  x3du2  A  dt 

da4  =  —  a2  A  a3  +  (x3  +  xi«2  —  1)<*3  A  dt  —  (x2u2  +  «i)a2  A  dt 

and  we  note  that  is  not  integrable,  thus  the  system  is  not  linearizable  by 

static  state  feedback. 

Now  consider  a  prolongation  by  differentiation  of  /  on  R10, 

/={J,/?\/?2,/?3} 
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with  the  additional  one-forms 


P1  =  du2  —  vdt 

ft2  =  dv  —  wdt 

P3  =  dw  —  z  dt 

and  additional  coordinates  v,wyz.  A  complement  to  /  is  {dui,dz,dt}.  The  derived 
flag  of  the  extended  system  has  the  form 

/=  u,\p\p\p*} 

/(1)=  »\P\P2} 

7<2>  =  U,4,  pi} 

/<3>  =  {u4} 

J<4>  =  {0} 

where  the  one-forms  adapted  to  the  derived  flag  are 

w1  =  dxi  —  u2dx2  +  (u22xi  —  x2)d< 

(J2  =  dx2  —  (u2Xj  +  x3)d< 

u>3  =  dx3  —  («i  -f  u2x2)dt 
u>4  —  dx4  —  iiidt 

The  structure  equations  are 

dw1  =  (—1  +  u23  +  »)  w1  A  dt  mod  1^ 

dw2  =  — u>3  A  dt  mod  1^ 

du>3  =  — dui  A  dt  mod  I 

dw4  =  -/31  A  dt 

dp1  =  -02Adt 

dp2  =  -03Ad< 

d03  =  -dzAdt 

from  which  it  is  easily  seen  that  each  {H'\dt}  is  integrable,  hence  the  prolonged 
system  I  can  be  feedback  linearized. 

We  note  that  although  this  example  can  be  linearized  by  dynamic  extension,  it  is 
shown  in  [4]  that  it  does  not  satisfy  the  sufficient  conditions  given  in  that  paper. 

Remark  2.  Dynamic  State  Feedback. 

We  have  shown  that  a  prolongation  by  differentiation  of  a  control  system  corresponds 
to  a  dynamic  extension.  A  general  dynamic  state  feedback  corresponds  to  adding 
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some  states  to  the  system  and  putting  feedback  around  them,  such  as 

x  =  f(x,u) 
z  =  g(x,z,v) 
u  =  P(x,z,v) 

This  general  form  does  not  correspond  to  a  Caxtan  prolongation,  since  there  may 
not  be  a  one-one  correspondence  between  trajectories  of  the  extended  system  and 
trajectories  of  the  original  system.  This  is  especially  obvious  if  the  added  z  states 
have  their  own  dynamics,  independent  of  x,  such  as  in  the  example  due  to  [12] 

x  =  /(x,u) 

=  *2 

u  =  P{z)v. 

where  harmonics  are  added  in  the  dynamic  state  feedback.  This  example  corresponds 
to  a  prolongation  of  the  original  system,  but  not  a  Cartan  prolongation,  since  there 
are  many  possible  trajectories  in  (x,  z)  space  for  every  trajectory  in  x,  depending  on 
the  initial  conditions  of  the  z  coordinates. 

7.  Conclusions 

In  this  paper,  we  presented  the  method  of  exterior  differential  systems  for  analyzing 
nonlinear  systems.  We  have  given  necessary  and  sufficient  conditions  for  converting 
Pfaffian  systems  to  Goursat  normal  form,  and  we  contributed  sufficient  conditions  for 
converting  systems  to  Goursat  form  using  prolongations.  In  addition,  we  showed  how 
the  techniques  that  we  described  for  general  Pfaffian  systems  could  be  specialized  to 
control  systems,  and  the  conditions  for  exactly  linearizing  systems  could  be  restated 
in  the  language  of  forms.  Since  dynamic  extension  is  the  dual  of  prolongation  by 
differentiation,  our  theorem  for  converting  Pfaffian  systems  to  Goursat  form  using 
prolongation  could  be  specialized  to  give  conditions  for  converting  control  systems  to 
Brunovsky  form  using  dynamic  extension  and  nonlinear  feedback.  We  showed  that 
these  conditions  are  closer  to  necessary  and  sufficient  than  those  which  exist  in  the 
literature. 

Future  directions  of  research  include  investigating  other  types  of  prolongation  be¬ 
sides  prolongation  by  differentiation  as  well  as  algorithms  for  systematically  prolong¬ 
ing  Pfaffian  systems  to  achieve  equivalence  to  Goursat  forms. 
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